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^ ■ ABSTRACT 
G ' 

^ ' We analyze a sample of ^^2 600 MlVS/Spitzer 24 /zm sources brighter than ~80/iJy and located 

1-^ ' in the Chandra Deep Field South to characterize the evolution of the comoving infrared (IR) energy 

, density of the Universe up to z ^ 1. Using published ancillary optical data we first obtain a nearly 

■ complete redshift determination for the 24 /im objects associated with i?<24mag counterparts at 
\ z <1. These sources represent ~ 55-60% of the total MIPS 24 /im population with /24pm ^ 80 /iJy, the 

■ rest of the sample likely lying at higher redshifts. We then determine an estimate of their total IR 
^ ' luminosities using various libraries of IR spectral energy distributions. We find that the 24 /xm pop- 

C\J I ulation at 0.5<z<l is dominated by "Luminous Infrared Galaxies" (i.e., lO"'^ Lq < Lir < 10"'^^ Lq), 

■ the counterparts of which appear to be also luminous at optical wavelengths and tend to be more 
\l ' massive than the majority of optically-selected galaxies. A significant number of fainter sources 

(5x 10^*^ L0 < Lir < 10^^ Lq) are also detected at similar distances. We finally derive 15 /im and total 
IR luminosity functions (LFs) up to z~l. In agreement with the previous results from ISO and 
SCUBA and as expected from the MIPS source number counts, we find very strong evolution of the 
contribution of the IR-selected population with lookback time. Pure evolution in density is firmly 
' excluded by the data, but we find considerable degeneracy between strict evolution in luminosity and 

^p^, a combination of increases in both density and luminosity (ijpj oc (l+z)"^'^-"-^ , oc {l+z)^'''-°-<i). A 

Q ' significant steepening of the faint end slope of the IR luminosity function is also unlikely, as it would 

| H I overproduce the faint 24 /im source number counts. Our results imply that the comoving IR energy 

c/3 . density of the Universe evolves as (l-t-z)^'^^" '* up to z ^ 1 and that galaxies luminous in the infrared 

■ (i.e., LjR > 10^^ L0) are responsible for 70±15% of this energy density at z~ 1. Taking into account 
^ I the contribution of the UV luminosity evolving as (l-|-z)""^-^, we infer that these IR-luminous sources 

. dominate the star-forming activity beyond 2~0.7. The uncertainties affecting these conclusions are 

' largely dominated by the errors in the fc-corrections used to convert 24 fim fluxes into luminosities. 

^ . Subject headings: galaxies: high-redshift — infrared: galaxies — cosmology: observations 



o 
in 
o 



1. INTRODUCTION 

The succ essful launch of the Spitzer Space Telescope 
ijWerner et al. 2004) recently opened a new exciting win- 
dow on the deep infrared (IR) Universe. Spitzer oper- 
ates between 3.6 and 160 /im with unprecedented sensi- 
tivity and better spatial resolution compared to previ- 
ous infrared satellites (e.g., IRAS, ISO). Directly prob- 
ing the dust emission or the redshiftcd signature of 
distant stellar populations, its first extragalactic sur- 
veys have already unveiled a huge numb er of faint 
and hi gh redshift source s (see f or instance IFazio et al.l 
2p04'; 'Eisenhardt et al.' '2004'; 'Paoovich et al.' '2004 
Mar le an et al. 2004; Lonsdale et al._, 2004; Char v et aQ 
20041 lYan et a l. 2004; Dole et al.ll2004al in the "Spitzer 
Special Edition - volume 154" of the ApJ Supplement). 
Spitzer therefore provides new opportunities to deter- 

^ Based on observations made with Spitzer, operated by the Jet 
Propulsion Laboratory under NASA contract 1407. 



mine the IR properties of galaxies in the general context 
of cosmic evolution. 

High redshift sources detected by Spitzer in the 
mid- and far- infrared wavelength range (i.e., 8 /im 
< A < 1000 /im) are characterized by intrinsically very 
high lumino sities llEg-ami e t al.ll2004 | iFrayer et al.ll2?)(il: 
llvison et all l2004t iLe Floc'h et all I2004D . They ap- 
pear as the distant analogs of the local Luminous and 
Ultra-Luminous InfraRed Galaxies (respectively LIRGs: 
10" Lq < LiR =L[8 - lOOO/iw] < 1012 L^, and ULI RGs: 
Ltr > 10^2 Lq, see the review by [Sanders & M irabel 
11996ft . Such infrared-luminous sources^ emit the bulk 
of their energy as dust-reprocessed thermal IR emis- 
sion powered by embedded star formation or by ac- 
creted material surrounding supermassive black holes. 
They were first discovered in the nearby Universe with 

We adopt in this paper the more general expression of 
"infrared-luminous galaxies" to denote sources characterized by 
LiR > 10" Lq. 



2 



E.Le Floc'h et al. 



ground-based observations ijRieke fc Low! 1197^. After 
being systematically catalogued by IRAS ijSoifer et all 
Il987|) . they were found to be locally very rare and to 
only account for ~5% of the t otal infrared energy emitted 
by galaxies at low redshift (Soifcr fc Neugebaueilll991t 
iKim fc Sanders,,,1998..) . Nevertheless, there is clear evi- 
dence that they were significantly more numerous ear- 
lier in cosmic history. In the past few years, deep 
observations performed in the infrared by ISO and in 
the submillimeter by the SCUBA camera have revealed 



strong 
time 


' evolution of these luminous sources with lookback 
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iDole et al. 
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apparent in the population of ra dio sources at /zJy 
flux levels fe.g.. ICowie et al"! 12004(1 . Characterized by 
a high redshift space density several orders of magni- 
tude larger than predicted by non-evolving models, in- 
frared luminous galaxies contribute a significant fraction 
of the distant starbursting activity and play a crucial 
role in the formation of massive spheroidals throughout 
the cosmic ages (e. g., Fl orcs ct al^ 1999; Gispcrt ct aL 
2000t IrVan ccschin i etall ]200lt Icharv fc Elbad |200IF 
Blain et aL.2002; .Chapman et al.ll2n08aD . 

This strong evolution of infrared-s elected sources has 
also been clearly seen by Spitzer fCharv ct al.' "2001, 
[Fapovich c t al. 2004 hereafter P04; Marlcau ct al. 2004; 
iDole et ani2004a}l . One of the most interesting results 
of the Spitzer deep surveys is the behavior of the dif- 
ferential number counts at 24 /itm. These counts turn 
over at fluxes lower than had been expected based on 
"pre-launch" models. Lag achc, e t al. (2004^ suggest that 
these counts reveal even more luminous galaxies z> 1.5 
than expected, though one could argue that they can 
also originate from a steeper faint-end slope of the in- 
frared luminosity function at more modest redshifts. To 
better understand the nature of the sources responsible 
for this turn-over and also more generally the role of in- 
frared galaxies in cosmic evolution, we examine in this 
paper the evolution of the comoving IR energy density 
with redshift up to z ~ 1. This study is based on a sam- 
ple of 24 /im-selected Spitzer sources within the Chandra 
Deep Field South and characterized by redshifts taken 
from the literature. A companion publication by Bell et 
al. (2005) explores in more detail the role played by these 
sources in the decline of the star formation history since 
2~0.7. 

The paper is organized as follows. In Sect. 2 we de- 
scribe the infrared and optical data used in this study, 
while Sect. 3 outlines the results of the cross-correlations 
that we performed among catalogs to determine the red- 
shifts of infrared-selected sources. In Sect. 4 we study the 
contribution of these infrared galaxies to the total counts 
at 2Afim for various redshift limits and present a com- 
parison with the predictions from various recent models 
of IR galaxy evolution. Using different libraries of tem- 
plates published in the literature, we further derive in 
Sect. 5 an estimate of the total infrared luminosities of 
these sources based on our mid-infrared observations. In 
Sect. 6 we analyze a few properties of their optical coun- 
terparts and we finally explore in Sect. 7 how the evolu- 
tion of the infrared luminosity function at < z < 1 in 
the CDFS can be constrained using our data. Interpreta- 
tions are discussed in Sect. 8, and we give our conclusions 
in Sect. 9. Throughout this work, we assume a ACDM 



cosmology with Hp = 70 km s~^ Mpc~^, rim = 0.3 and 
r^A = 0.7 (jSoergel et alJl2003D . Unless exphcitly stated, 
magnitudes are quoted within the Vega system. We also 
adopt the universal Initial Mass Function from lSalpeteil 
1IT95I . 

2. THE DATA 

2.1. 24 fJ-m imaging 

The region of the Chandra Deep Field South (here- 
after CDFS, a = 3''32™00^ S = -27°35WlJ2000}_was 
observed with the MIPS instrument l|Rieke et al J 120041) 
on-board the Spitzer Space Telescope in January 2004 
as part of the MIPS Guaranteed Time Observing pro- 
gram. These observations were performed over a total 
field of -1.45x0.4= 0.6 deg2 with the so-called "Slow 
Scan" technique, a MIPS observing mode that allows the 
coverage of large sky areas with high efficiency. The de- 
tector at 24 /im uses a 2.45" pixel size array of 128x128 
elements and the image at this wavelength is character- 
ized by a Point Spread Function (PSF) with a Full- Width 
at Half Maximum (FWHM) of -6". The effective inte- 
gration time per sky pixel was —1380 s. Data reduc- 
tion w as carried out using the MIPS Data Analysis Tool 
(DAT. lGOTdbn et al. 2005). The mosaic finally obtained 
has a pixel scale of half that of the physical detector 
pixel. A color version of this 24 /im map can be seen in 
the publication by Rigbv et al. (2004|. A sub-region of 
this field is also illustrated in Figure 1. 
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Fig. 1. — A sub-image of the CDFS 24 /im observations obtained 
with MIPS, centered at a = 3''32™28= and S = -27°48'27" 
(J2000). The fields of view respectively covered by COMBO-17 
(dashed line), VVDS (solid line) and GOODS (dotted line) are 
also indicated (see text for details) . The 24 /im sources consid- 
ered in this paper are located within the 775 arcmin^ of overlap 
between MIPS and COMBO-17. 

Since most sources are unresolved in our data, ex- 
traction and photometry were performed using the PSF 
fittin g technique of the DAOPHOT software ijStetsonI 
Il987t) . An empirical point spread function was con- 
structed from the brightest objects found in our mosaic, 
and it was subsequently fitted to all the sources detected 
in the map. Allowing for multiple-match fitting to deal 
with blended cases, we derived the flux density of each 
source from the scaled fitted PSF and finally applied a 
slight correction to account for the finite size of the mod- 
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eled point spread function. We also performed extensive 
simulations adding and recovering artificial sources in 
the data, which allowed us to derive an 80% complete- 
ness limit at ~83 /xJy. Contamination by false sources 
appears in our detection catalog at fluxes fainter than 
90 /LtJy. In the full sample of objects brighter than the 
80% completeness limit, we estimate such a contamina- 
tion to be around 0.6%. A more de tailed de s cripti on of 
our technique is presented by Papovi ch et al.l lj2004f) who 
also discuss the evolution of the completeness limit and 
the false source fraction as a function of the flux at 24 fim 
(see their figure 1). 

Distortion effects are efficiently removed by the DAT. 
The final astrometric uncertainty of the mosaics pro- 
duced by the pipeline therefore originates from the point- 
ing reconstruction. To quantify this uncertainty in our 
data, we cross-identified the brightest sources detected 
at 24 /xm with the Two Micro n All Sky Survey (2MASS) 
catalog l|Jarrett et al.l I200CI|) . A systematic offset of 
^0.6" was observed in the scan direction^ and subse- 
quently removed for our source catalog to match the 
2MASS coordinates. We estimate that the remaining 
scatter relative to 2MASS is better than 0.3" rms. 

2.2. Optical- spectroscopic redshifts 

To derive the redshifts of the MIPS 24 /xm-selected 
sources, we retrieved from the literature publicly avail- 
able catalogs of optical spectroscopic surveys in the 
CDFS such a s the "VIMOS VLT D eep Survey" (here- 
after VVDS, ILe Fevre et al] l2?ffll . the ES0/F0RS2 
survey performed by the GOODS Legacy team^ 
ijVanzella et 3.1.1200511 . and the foll ow-up of X-ray sources 
described bv lSzokolv et all l|200l 6 ^he VVDS consor- 
tium has released redshifts for 1599 sources located in an 
area of 21x21.6arcmin^ with an overall redshift mea- 
surement completeness of ^^88% down to Iab < 24. 
For each identification, a flag indicates the reliability 
level of the corresponding measurement; 1457 sources 
in this survey are classified with more than 75% con- 
fidence in the redshift determination. The GOODS cat- 
alog provides 234 redshifts obtained within a slightly 
smaller area of the CDFS (lOx 15 arcmin^). They are 
also tentatively classified into three categories depend- 
ing on their reliability, and 150 sources in this catalog are 
thought to have a se cure redshift measurement. Finally, 
iSzokolv et al.l l)2004|l present a spectroscopic follow-up of 
Chandra X-ray sources, with redshifts obtained for 168 
objects. Among those, 126 are considered to have unam- 
biguous identifications. 

We combined the catalogs of these three optical sur- 
veys to create a single list of 1 941 spectroscopic redshifts 
with their corresponding flags. Because of some overlap 
between the different observed regions, several redshifts 
were sometimes assigned to a given single source. In case 
of discrepant estimates, we kept the one flagged with the 
highest confidence in the redshift measurement. 

The MIPS image at 24 /im entirely covers the areas 

^ This offset seems to be due to systematics related to the po- 
sition of the MIPS scan mirror. It affects the headers of the data 
obtained before May 2004 {Spitzer Science Center, priv. commu- 
nication) . 

* http:/ /cencosw. oamp.fr/EN/index. en. html 

^ http:/ /www. eso.org/science/goods/spectroscopy/products. html 
® http://www.mpe.mpg.de/ mainieri/cdfs_pub/ 
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observed by these spectroscopic surveys. The VVDS and 
GOODS fields of view can be seen in Figure 1. 

2.3. Optical-photometric redshifts 

In addition to the spectroscopic redshifts previously 
described, we also made extensive use of the photo- 
metric redshifts from the COMBO-17 survey {"Clas- 
sifying Objects by Mediu m-Band Observations in 17 
filters", IWolf et abl l200l . COMBO-17 observed a 
31.5x30arcmin^ region of the CDFS (see Fig. 1) through 
a set of 5 broad-band and 12 narrow-band filters, which 
allowed the determination of accurate spectral energy 
distributions (SEDs) and absolute magnitude s for several 
thousand optically-selected sources (see also iWolf et alJ 
'2003*). Each SED was analyzed using a library of repre- 
sentative templates for various spectral types, and a red- 
shift probability distribution was subsequently derived 
for each source. For a total of 24 217 objects, these 
distributions were successfully fitted with Gaussian-like 
functions, the mean and the variance of which led to 
well-constrained estimates of photometric redshifts with 
corresponding uncertainties. 

These redshifts are accurate to 1% in (5^/(1 -I- z) for 
galaxies with R < 21 mag, and they are mostly reli- 
able (i.e., Sz/{1 + z) < 10%) for all sources at z<1.2 
and brighter than i?~24mag (11422 objects). The 
yield of high weight redshifts drops steeply for 2:>1.2 
or i?>24mag, so we did not consider such faint or dis- 
tant objects. Another redshift estimate is also provided 
for nearly every detection of the survey (62 366 sources) 
based on the peak of the computed probability distribu- 
tion. B ecause it is less r eliable, it should be used with 
caution l)Wolf et al.l2004tl and was therefore also not con- 
sidered in this work. About '^85% of the region covered 
by COMBO-17 overlaps with our MIPS 24 //m mosaic 
and the spectroscopic optical surveys. This overlapping 
area has a total field of view of 775arcmin^ (see Fig. 1). 

3. OPTICALLY-SELECTED REDSHIFTS OF MIPS 
SOURCES 

3.1. Cross- correlation between the infrared and the 
optical data 

From the full MIPS catalog we first selected the 24 fim 
sources located in the common area covered by Spitzer 
and the aforementioned optical redshift surveys (i.e., 
775arcmin^). In this region of overlap we detected 
5 589 objects at 2A^m, with 3 616 of those having a flux 
greater than the 80% completeness limit of 83 /iJy. 

We cross-correlated this 24 iim sub-sample with the op- 
tical data using a tolerance radius of 2" for matching 
sources. This choice was first motivated by the rather 
large FWHM of the MIPS 24 ^m PSF (~6") compared 
to that typically seen in optical images. As already ob- 
served in local interacting systems (e.e..\Le Floc'h et alJ 
[200l ICha,rma,nda,ris et a].ll2004t IGa,l]a,is et a,].ll2004^. it 
also accounts for the physical shift that could be present 
between the location of the infrared emission and the 
brightest optical component of distant mergers (2" cor- 
responds to a linear projection of ~15kpc on the sky 
at z~l). It is yet reasonably small given the position 
accuracy of the source centroids at 24 /im (better than 
0.5" rms). Taking a larger value also increases the risk 
of associating MIPS sources with wrong optical counter- 
parts in the case of multiple matches. When a double 
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match was found, we selected the closest object. We ig- 
nored the cases where three or more optical sources could 
be associated with a given 24 /im detection. 

We first correlated the 24 /im source catalog with the 
list of spectroscopic redshifts. A total of 543 matches was 
found in this cross- identification, of which 465 objects are 
flagged to have a high-confidence redshift measurement. 
The fraction of multiple matches was only 1%. At this 
stage, we only kept the identifications with a secure red- 
shift determination. Those 465 sources represent only 8% 
of the infrared-selected sub-sample, which emphasizes 
the critical need for using photometric redshifts. We thus 
cross-correlated the rest of the data with the catalog of 
COMBO-17. We found 2 170 MIPS sources (1 987 single 
and 183 double matches) identified with a photometric 
redshift below z = 1.2 and an optical counterpart brighter 
than i? = 24mag. For 9 objects selected at 24 /im (less 
than 0.5% of the sample) three possible matches were 
found within 2" around the MIPS source. They were 
not further considered. In total, we assembled a cata- 
log of 2 635 MIPS sources (of which 1 962 are brighter 
than the 80% completeness limit of the 24 /im survey) 
associated either with a reliable spectroscopic redshift or 
a clearly-constrained photometric redshift. Virtually all 
are at z < 1.2, since the yield of values at higher redshift 
is very low with both spectroscopy (due to the "redshift 
desert") and COMBO-17. We also identified another 
set of 1 681 MIPS sources with an optical counterpart 
in the COMBO-17 catalog but without any reliable red- 
shift (271 of them have /24pm > 83 /iJy and are brighter 
than i? = 24mag). 

3.2. Redshift uncertainties 

Since only a very small fraction of the MIPS sources 
have been identified with secure spectroscopic redshifts, 
it is worth looking at the typical uncertainties of the 
other photometric redshifts characterizing the 24 /im- 
selected sources. The accuracy of the COMBO-17 clas- 
sification decreases for sources fainter than i? ^ 22 mag. 
This may have non-negligible effects when estimating 
e.g., source densities as a function of lookback time, espe- 
cially when the uncertainties become comparable to the 
redshift bins in which galaxy properties are averaged. 
Given that the 24 /im sample is by definition selected 
through the emission by warm dust, one may further- 
more question whether the implied extinction at optical 
wavelengths can lead to a more significant redshift mis- 
classification in the specific case of the most luminous 
(i.e., dust-obscured) MIPS sources. 

In Figures 2a & 2b we compare high-confidence spec- 
troscopic redshifts of optically-selected field galaxies and 
MIPS 24 /im detections with their photometric redshifts 
estimated by COMBO-17. We see that the photo- 
metric redshift errors are small (i.e., |zspcc-^phot| ^ 0.1) 
and they are not statistically larger in the case of the 
MIPS sources. The latter can be explained as follows. 
Mid-infrared space-borne and ground-based observations 
of local LIRGs/ULIRGs reveal that the dust respon- 
sible for the bulk of the IR luminosity of those ob- 
jects originates from very c ompa,ct reg ions l)Soifer et al.l 
I2OOO, 2001; Charmandaris eTal] l200^ . The effect of 
extinction in these dusty systems is therefore very lo- 
calized and is usually not apparent in the global spec- 
tral energy distribution of their optical counterparts 



(jSanders fc MirabeMl996fl . A similar situation likely per- 
tains at higher redshifts up to at least z^l. Atz'^0.7 
for instance, ISOCAM 15 /xm and MIPS 24 /im-selected 
galaxies are indeed associated with luminous optical 
sources characterized b y a wide range of optical colors 
and r norphologies ( Flores et al]ll999t iRigopoulou et alJ 
120021 iFranceschini et al.ll2005ir Bell et al. 2005; see also 
Sect. 6). Only a careful approach based on medium- 
resolution spectroscopy can distinguish the m from the 
optic ally-selected galaxy population (e.g., F lores et alJ 
|2004|) . It is therefore unlikely to encounter any signif- 
icant increase of redshift mis-classification as a function 
of infrared luminosity (at least up to z~ 1). 

Finally, Fig. 2c shows the photometric redshift uncer- 
tainties provided by the COMBO-17 catalog as a func- 
tion of the observed i?-band magnitude for the MIPS 
sources considered in this paper. These errors were de- 
termined as the square root of the variance characteriz- 
ing the redshift probability distribution of each object. 
In agreement with the comparison we made from Fig. 2a 
using the spectroscopic sub-sample, we see that sources 
brighter than mag have a very accurate redshift 

estimate (i.e., Sz < 0.04) and most of the sample (88%) is 
characterized by a redshift accuracy better than (5^ ~ 0.1. 
The average uncertainty is only <5z>— 0.07 with a dis- 
persion (72 =0.05 for sources with 22mag<i? < 24 mag. 
It rises to <Sz>=0.09 (with a similar dispersion) if we 
only consider objects fainter than i? ~ 23 mag. Looking 
at the evolution of the MIPS sources averaged per red- 
shift bin of Az^O.2 should lead therefore to rather ro- 
bust results. 

3.3. Completeness of the sample 

The biggest concern affecting our results is the incom- 
pleteness of the redshift determination for the most dis- 
tant sources considered in this paper. As previously men- 
tioned, 3 616 objects were detected by MIPS above the 
80% completeness limit in the region overlapping with 
the optical surveys of the CDFS, while our final catalog 
only contains 1 962 sources with f 24,^1111 > 83 //Jy and a 
redshift identification (i.e., ^55%). To better assess this 
effect, we plot in Figure 3. a the fraction of MIPS sources 
characterized either by a spectroscopic or a photometric 
redshift as a function of the observed 24 /im flux (solid 
line) as well as the fraction of their identification with an 
optical counterpart for various magnitude limits (shaded 
regions). At the highest fluxes (i.e., /24^m <L IniJy), 
the incompleteness of the sample is mostly due to a few 
objects detected close to very bright stars blooming the 
optical image and where a reliable identification of the 
counterpart cannot be obtained. These cases should not 
introduce any bias in our results. Below 1 mJy, the drop 
in the redshift determination completeness corresponds 
to sources fainter than i?^ 23.5 mag (see also Fig. 3.b). 
Above the 80% completeness limit of our MIPS data (i.e., 
83/iJy), we see for instance that a significant fraction 
(^ 25%) of the MIPS objects with no redshift are associ- 
ated with 24 mag <R< 25.5 mag sources, while ~ 20% of 
them have counterparts fainter than i?> 25.5 mag. This 
also indicates that a complete study of the faint MIPS 
population will likely require an extensive use of photo- 
metric redshift techniques. 

The consequence of this limitation can be seen by com- 
paring the redshift determination completeness as a tunc- 
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Fig. 2. — a) Comparison between the VVDS spectroscopic and COMBO-17 photometric rcdshifts as a function of the i?,-band magnitude, 
for 1142 optically-selected galaxies (dots) and 308 MIPS 24 ^m sources (open squares), b) Histogram of the absolute values Izspcc-Zphot 1 1 
with the optically-selected and 24 ^m sources respectively indicated by the solid and dashed lines, c) Photometric redshift uncertainties 
provided by COMBO-17 for the 24 fini sources considered in this paper. Selection cut is set to < 24 mag and z < 1.2. In agreement with 
comparisons from panels a) and b), the accuracy is better than Sz ~0.1 for 88% of the sample. 



tion of the i?-band magnitude with the distribution of 
the i?-band magnitudes of the MIPS sourees for dif- 
ferent redshift limits. This is shown in Fig. 3.b, while 
Fig. 3.C represents the absolute i?-band magnitudes of 
the MIPS sources as a function of redsfhit and derived 
from COMBO-17. We see that the identification should 
be nearly complete up to but we start missing 

24 yum objects located at higher redshift and associated 
with optical counterparts fainter than JS~ 23.5-24 mag. 
These mis-identifications likely induce a biais against the 
faintest sources detected at 24 /um. In Sect. 6 we will 



establish a correlation between the optical and the in- 
frared luminosities of galaxies in our sample, and we will 
use such relation in Sect. 7 to quantify how this bias af- 
fects the estimates of luminosity functions (see also the 
Appendix for further details characterizing these faint 
24 /im sources). However, given the importance of the 
cosmic variance arising from the variations of large scale 
structures in the CDFS (see Sect. 3.4), and taking into 
account the errors in the ^-corrections used to derive 
IR luminosities (see Sect. 5), we infer that the possible 
missing redshift determinations should not dominate the 



6 



E.Le Floc'h et al. 



absolute uncertainty in the source density estimate up to 
z-1. 
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Other arguments also suggest that the bulk of the 
MIPS sources for which we could not identify any re- 
liable redshift should be indeed located at z > 1 and will 
not affect this study. First we will show in the fol- 
lowing sections that infrared luminous galaxies are as- 
sociated with optical counterparts at the bright end of 
the luminosity function in the visible. Luminous opti- 
cal sources with L > Lt^ but fainter than i?~ 23.5-24 mag 
should lie at z > 1 (Bcnitcz 2000). Moreover, preliminary 
results from Spitzer reveal that a significant fraction 
of the MIPS sources are located at such large redshifts 
iCharv et all l200l lEgarni et all l200l ILe Floc'h etUT] 
12004 iLaeache et al.ll2004j) . Finallv. although they are 
not fully reliable, the COMBO-17 redshift estimates 
based on the peak of the redshift probability distribu- 
tions indicate that the majority of the sources that we 
have not identified should be at z'>l. 

3.4. Redshift distribution 

Figure 4 shows the distribution of the redshifts derived 
for our sample of infrared sources (bottom panel, solid 
line). It is compared to the redshift distributions of the 
VVDS (dotted line) and the COMBO-17 objects with 
i? < 24 (dashed line). As previously discussed, the 
redshift identification of the MIPS sources is complete 
only up to z ^ 0.8, and the distribution beyond this limit 
should thus be considered as a lower limit estimate. 
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Fig. 3. — a) Fraction of MIPS sources identified with a spec- 
troscopic or a photometric redshift as a function of the observed 
flux at 24 /^m (solid fine). We also indicate the fraction of MIPS 
sources with an optical counterpart detected in the i?-band for 
various magnitude limits (shaded regions). The vertical dotted 
line corresponds to the 80% completeness limit of the 24 fira ob- 
servations, b) Fraction of 24 /im source optical counterparts iden- 
tified with a redshift, as a function of the i?-band magnitude (solid 
line). Shaded regions show the i?-band magnitude histograms (all 
scaled with an arbitrary constant factor) of the MIPS sources up 
to various redshift limits. The drop in the redshift identification is 
clearly apparent at H ~ 23.5-24 mag. The sample should be com- 
plete up to 2 ~ 0.8. c) Absolute ij-band magnitudes of the MIPS 
sources as a function of redshift. 



Fig. 4. — Bottom: The redshift distribution of 24 /im- 
selected sources over a 775 arcmin^ region within the CDFS 
(solid line), compared to the distributions of photometric red- 
shifts for R< 24 mag sources in the COMBO-17 catalog (dashed- 
line, scaled down by a factor of 2.5) and spectroscopic redshifts 
for Ia_b < 24 mag sources in the VVDS survey (dotted line). The 
contribution of MIPS sources with spectroscopic redshifts is also 
indicated (dash-dotted line). Note the prominent overdensity at 
2^0.65, clearly seen in all the distributions. Source densities at 
z > 0.8 should be considered as lower limit estimates due to the 
incompleteness of the redshift identification. Top: The fractions 
of R < 24 mag sources detected at 24 fim, estimated as a function 
of redshift in the whole population of COMBO-17 (dashed line) 
and restricted to the objects with Mb <— 16 (solid line). 
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We note that the three distributions (i.e., MIPS, 
VVDS, COMBO-17) are more or less similar from z~0 
to z~l. Up to this redshift limit, MIPS is indeed sensi- 
tive to sources luminous in the infrared but also to more 
normal galaxies (see Sect. 5) such as those detected at op- 
tical wavelengths and driving the VVDS and COMBO-17 
redshift distributions. Furthermore, we note a significant 
number of objects around z~0.65, apparent in all the 
distributions. This originates from an overd ensity char- 
acteristic of the CDFS near this redshift fe.g.. lWolf et al.l 
I2OO4') and it is likely rel ated to cosmic var iance and large- 
scale structure effects ()Somerville et al. 200^. It may 
constitute a non negligible source of uncertainty in our 
density estimates. A comparison between the integrated 
i?-band luminosity densities as a function of redshift from 
galaxies of the blue sequence (Bell, private communica- 
tion, see Willmer et al., in prep., for the evolution of the 
blue sequence luminosity function) in the CD FS and av- 
erage d over the three fields of COMBO-17 l|Wolf et all 
1200,'^ shows that the CDFS is most often underdense 
by ~50% at 0< z <1, except indeed at z-0.65±0.1 
where the overdensity reaches ~20%. Given the simi- 
larity between the redshift distributions of the optically 
and infrared selected sources, and since most of MIPS 
sources at z ^^0.7 appear to be associated with large spi- 
rals dominating the i3-band emission (Bell et al. 2005, 
see also Fig. 10b), this cosmic variance characterized from 
the i?-band luminosity density is likely affecting also the 
population of galaxies detected at 24 /xm. Therefore the 
apparent peak a,t 0.65 in the MIPS redshift distri- 
bution is probably only related to the structure of the 
CDFS and it should not be a characteristic of the gen- 
eral population of infrared-selected sources. 

In spite of the similarities that we observe between 
the redshift distributions in the visible and the infrared, 
a more critical look reveals that the fraction of optical 
sources brighter than Mb = -16 and detected at 24 /im 
decreases from z = to z~0.4 and then increases with 
redshift up to Z'^l (see Fig. 4., top panel). Interpreting 
this trend is not straightforward, but it could reflect the 
difference in the fc-correction effects observed at optical 
and infrared wavelengths. As we will see in Sect. 5 (see 
Figs. 7 & 9) the effective sensitivity of the deep 24 /zm ob- 
servations drops rapidly from z = to z ^0.5, which ex- 
plains the sharp decrease of relative 24 /im detections in 
this redshift range. Then it remains roughly constant at 
0-5<z<0.9 because of strong emission features present 
at 11.3 and 12.7 /xm in galaxy spectra progressively enter- 
ing the 24 /j,m filter at these redshifts. Since the sensitiv- 
ity of the optical data constantly drops with distance, an 
increase in the relative number of 24 /im sources detected 
beyond z ~ 0.5 in the visible might thus result. Note that 
it may also reflect the stronger evolution of IR sources 
with lookback time, since IR-luminous p hases in galax- 
ies were more common in the past fe.g.. iHammer et al.l 
[2005). The lack of apparent decline moreover suggests 
that the intrinsic peak of the MIPS population redshift 
distribution must lie at least beyond 2;~0.8. 

Finally, another overdensity lies at 2:^0.15. It repre- 
sents a rather small fraction of our detections at 24 /im 
(i.e., 4%). Since we will only consider sources at z > 0.3 
when building the luminosity functions in Sect. 7, it 
will not affect our final results on the evolution of the 
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IR galaxy population. 

4. EXPLORING THE ORIGIN OF THE BREAK IN THE 
MIPS 24 fim NUMBER COUNTS 

4.1. Cumulative differential counts as a function of 
redshift 

Using our matched catalog of 24 /xm-selected objects 
with redshifts, we now explore for various redshift slices 
the contribution of these MIPS sources to the total differ- 
ential number counts. Our goal is a better understanding 
of the origi n of the bump detected at 0.2-0.4 mJy (P04, 
iMarleauet a l. 2004). Figure 5 illustrates these cumula- 
tive source counts for redshift limits of 0.3, 0.5, 0.8, 1.0 
and 1.2 along with the global contribution of the MIPS 
sources identified with an optical counterpart in the full 
COMBO-17 catalog (i.e., i?<25.5mag). The distribu- 
tions for redshift limits z = l and z = 1.2 are likely un- 
derestimated due to the incompleteness of the sample as 
showed in Sect. 3.3. However, uncertainties due to cosmic 
variance likely dominate here (~ 0.10-0.15 dex based on 
the ratio between the i?-band luminosity densities of the 
CDFS and averaged within the other fields of COMBO- 
17). For comparison, we also represent the total 24 /im 
number counts de rived from our s ub-sa mple as well as 
those obtained bv lPapovich et all l)2004D . who averaged 
the 24 /im source population over ~ 10 deg^ in several 
cosmological fields. For consistency with our data, these 
counts from P04 were not corrected for incompleteness 
of the MIPS detections at faint fluxes. 




0.1 1.0 

24/im flux (mJy) 



Fig. 5. — Cumulative differential number counts observed at 
24 /im as a function of redshift (shaded regions). Distributions 
at z > 0.8 should be considered as lower limit estimates due to 
the incompleteness of the redshift identification. The striped area 
represents the contribution of all 24 fira sources identified with an 
optical counterpart down to ij<25.5mag. The total differential 
counts derived from our sub-sample are illustrated by the b lack 
solid line, while those determined by IPapovich et al.( <2004l) are 
represented by the dashed line and the vertical error bars. Counts 
are normalized to the Euclidian slope. 

We note that the 24 /im number counts derived in the 
CDFS show a roughly good agreement with those ob- 
tained by P04 below /24pm ~ 0.2 mJy. However, the den- 
sity of 24 /im CDFS sources at higher fluxes appears to 
be slightly lower (^20%), and consequently the peak of 
the differential number counts seems to occur at fainter 
levels. As previously seen in Sect. 3.4., such variations 
are easily understood in terms of cosmic variance and the 
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smaller number of infrared-bright sources in the CDFS is 
likely related to the underdensity observed in the _B-band 
at z<l. 

At critical fluxes f24fj,m ^ 0.2-0.4 mJy where the 24 fim 
differential counts present a turn-over, sources lying at 
z <1 contribute ^ 55% to the whole sample. Even 
though this fraction could be slightly higher due to 
the possible incompleteness of our redshift identification 
at 0.8<z<l, a significant part of the MIPS popula- 
tion should therefore be located at higher redshifts (i.e., 
z>l). As already derived from Fig. 3, it can also be 
noted that ~ 20% of the 24 /xm sources have optical coun- 
terparts fainter than 25.5 mag. 

4.2. Comparison with model predictions 

In Figure 6 we compare the differential counts in four 
redshift slices between z = and 2 = 1.2 with the pre- 
dictions fr om the models of IR gal axy evolut io n pro - 
nosed bv ILa.gacheet a1.l (^QM), Cha,rv et al.l ^2004) . 
Grunpioni et al.l l)200,'i[ see also Pozzi et al. 2004) and 
PearsonI l|2005j) . These models are all tied to the local IR 
galaxy population but differ in (i) the description of its 
components and their global properties (i.e., SEDs, lumi- 
nosity functions per object type) and (ii) the assumptions 
and parameterization used to infer the backward evolu- 
tion of IR sources. They provide a reasonably good fit to 
the total 24 /im number counts. They also reproduce a 
variety of other observables such as the IR background, 
the counts and/or the redshift distributions of the ISO 
and SCUBA galaxy populations. 
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Fig. 6. — Differential 24 /xm number counts produced by the 
MIPS sources in four redshift shces between z = Q and 2 = 1.2 as 
indicated within each panel (shaded reg ions). They are com- 
pared to the prediction s bv ILagache et a l. (2004, dashed lines), 
IGruDDioni et al.l 12003. dotted lines). ICharv et al. (2004, dash- 
dotted lines) and IPearson ( 200f| triple dot-dashed line) for simi- 
lar redshift ranges. The total counts derived from our sample are 
shown in each panel by the black solid line. The vertical bar in the 
upper panels represents the typical uncertainty due to the effect of 
large scale structures. The counts at 0.8 <z< 1.0 and 1.0 < 2 < 1.2 
(lower panels) only show a lower limit given the incompleteness 
of the redshift identification in our sample at z > 0.8. 

Looking at their predictions for the 24 /xm num- 
ber counts within the redshift slic es considered in thi s 
Fig. 6, we see that the models from ICharv et all l)2004j) . 



iGruppioni et alJ l)2005j) an d|P earso ^ ll200l may underes- 
timate the contribution of faint sources at < z < 0.8 but 
predict too many bright objects at z > 0.5. ILagache et alJ 
|2004), on the other hand, closely follow the observations 
at 2;<0.5 but might globally underestimate the source 
density at 0.5 < 2 < 0.8. All these models finally seem to 
overpredict the contribution of bright galaxies at z > 0.8. 
However we also note that these comparisons could be 
undermined given not only cosmic variance effects and 
the global underdensity of the CDFS between z = and 
z ~ 1 but also the incompleteness of our redshift identi- 
fication at z > 0.8. 

Comparisons with other models in the literature such 
as those published before the launch of Spitzer lead to 
larger discrepancies. These models tend to predict the 
characteristic turn-over of the differential 24 /zm number 
counts at much higher fluxes than where it has been ob - 
served (see for instance Figure 3 of lPapovich et alJ2004j) . 

5. TOTAL INFRARED LUMINOSITIES 

5.1. Estimating the total IR luminosity of galaxies from 
their mid-IR emission 

The rest- frame mid-infrared (MIR) regime (i.e., 8 /im 
< A <40 /im) is considered to be a good tracer of the 
bolometric luminosity of galaxies. Using the 12 /im IRAS 
galaxy sample, ISpinogho et al. (1995) first showed that 
the relative dust content of galaxies balances their to- 
tal energetic output between the optical and the far- 
infrared, leaving a "pivot point" in the MIR where the 
specific luminosity linearly scales with the bolometric 
on e. More recently , a simi lar conclusion has been reached 
bv ICharv fc Elbaa l|200H) who found a tight correlation 
between the 15 /im and the total infrared luminosities 
for a sample of normal and luminous sources studied 
with ISOCAM. Based on these considerations and taking 
into account the general IR/submm color-color and/or 
luminosity-color correlations observed in the local Uni- 
verse, several authors have built libraries of luminosity- 
or color-dependent galaxy templates, which can then be 
used to estimate the total infrared luminosit y of galax- 
ies from their 24 /xm flux densities (Dale et al. 20011 
ICharv fc Elbadll2001t iDale fc Helo'ull2002; .Laeache et aU 
M03; Chanial 2003,). 

Figure 7 shows the total infrared luminosity that can 
be detected in our survey down to a 24 /im sensitivity 
limit of 0.08 mJy, as a function of redshift and computed 
with the aforementioned libraries. Provided these tem- 
plates are still representative of high redshift sources, we 
see that MIPS can easily detect normal starbursts (i.e., 
Ljr^IO^^Lq) up to Z'^1, and LIRGs up to z^l.5. 
Note that such limits only characterize a pure selection 
of sources at 24 /im, and the use of a priori informa- 
tion from detections at other wavelengths can allow the 
identification of MIPS sources with similar infrared lumi- 
nosities at even higher redshift (see e.g.,jLe Floc'h et aQ 
|2004() . The 80% completeness level r oughly co rresponds 

to a 5(7 detection in the MIPS data l|Dole et al 2004b.) . 

and the quoted 0.08 mJy limit is therefore rather conser- 
vative. 

Figure 7 also shows the dispersion between the pre- 
dictions (and th erefore the templ ates) of the various li- 
braries (see also iPapovich fc Bel]ll2002il . For a given ob- 
servable (e.g., infrared luminosity, color, ...) these li- 
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braries assign a unique SED that slightly varies from 
one set of templates to another. Such variations re- 
flect an intrinsic dispersion of the different correlations 
observed between the MIR/FIR/submm properties of 
galaxies th at is usually not taken into account (but see 
iChapmiaTet al 2QQ3b for the use of bi-variate luminos- 
ity functions addressing this issue). As an example, 
we compare in Figure 8 the observed relation between 
LjR and the monochromatic 15 /im luminosity Li5^ for 
a sample o f local galaxies studied with ISO and IRAS 
(jCharv fc Elbaz 2001i) with the predictions of the tem- 
plate libraries considered in Fig. 7. We see that the dis- 
persion can easily reach ~0.2dex. 
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Fig. 7. — Detection limits in terms of total infrared lumi- 
nosity, derived for two different MIPS 24 /.tm s ensitivities and as- 
suming the template libraries of Lagache et al. 12004, thiek solid 
line), Lagaehe et al. (2003, thin solid line), Dale et al. (2001, dot- 
ted linel. Chary fc Elbaa I.200K dashed line) and lChaniali t200a . 
dash-dotted line) as a function of redshift. The lower (respec- 
tively, upper) set of curves corresponds to a flux limit of 0.08 mjy 
(0.3 mJy) typical of the MIPS deep (shallow) surveys. 

To explain this dispersion, one may note that in the lo- 
cal Universe the mid-infrared SEDs of starburst sources 
with comparable bolometri c luminosities are s ubject to 
significant variations (e.g.. lArmus et aLil2004|) . In the 
rest-frame wavelength range probed by the MIPS 24 /i,m 
filter for high redshift galaxies, the MIR emission re- 
sults from the combination of prominent broad-band fea- 
tures mostly observed between 3 and 14 /im and usually 
denoted the Polycyclic Aromatic Hydrocarbons bands 
(PAHs), superimposed on a rising continuum of Very 
Small Grains (VSGs) stocha stically heated by th e young 
star radiation field (see e.g.. lLaurent et al]l2000l and ref- 
erences therein). PAHs present an amazingly universal 
SED signature as a global set of features wheii inte- 
grated over normal spira l galaxies ijRoussel et al.ll200lL 
but see lSmith et alJl2004l for smaller scale variations and 
a newly-discovered 17.1 /xm PAH in NGC7331). How- 
ever, they are not detected in low me tallicity sources 
l|Thuan et all 119991: Hfouck et al.ll2004al: Engelbracht et 
al. 2005, submitted). They are also believed to be de- 
stroyed within intense radiation fi elds as in the vicinity o f 
active galactic nuclei (AGN, e.g.. lLe Floc'h et al.ll200H) . 
Finally, the temperature of the VSGs as well as the sil- 

We define the monochromatic luminosity e&\ji, =v Y. Sv , with 
Sv the monochromatic flux of the galaxy expressed in WHz^^. 
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icate absorptions at 9.7 /im and 18 /im represent other 
important factors shaping the underlying continuum of 
the global mid- infrared SED of galaxies. All of these 
variations can therefore bring significant uncertainties in 
relating the MIR emission to the full infrared range of 
galaxy SEDs. 
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Fig. 8. — Ratio between the total IR and the monochromatic 
15 /im luminosities as a function of Lij^ predicted by the SED 
libraries considered in Fig. 7 (similar style coding). The grey area 
represent s the relation ( mean value with Icr boundaries) derived 
by Charv fc ElbaJ 12001]) for galaxies that were observed with ISO 
and IRAS. 

Furthermore, one may question whether these tem- 
plates derived from the properties of local galaxies are 
truUy representative of higher redshift sources. Infrared 
luminous galaxies detected by ISOCAM at a median 
redshift of ~ 0.7 are characteriz ed by roughly half-solar 
metallicites llLiang et al.l l2004f) , which could point to 
IR SEDs slightly different from those of local LIRGs. 
However, PAHs are still detected at this metallicity 
range in nearby starbursts (Engelbracht et al. 2005, 
submitted). Moreover, the 24/15/im flux ratios ob- 
served in 0~O.7 ISOCAM sources still a rgue for the 
prese nce of significant PAHs in their SEDs l|Elbaz et alJ 
l2005j) . and prominent PAH features were also recently 
obser ved with the Infra Red Spectrograph on Spitzer 
(IRS, 'Houck et al."20G4b") at even higher redshifts (i.e., 
1.7<z<2.8, Houck ct al. 2005). We also note that the 
MIR/radio relation expected from the MIR/FIR and 
FIR/radio correlations observed locally (lCo ndonlll992t 
iCharv fc E lbaz 200^ appears to be still v alid i n the more 
distant Universe (e.g.. lElbaz et alJl2002t [Appleton et alJ 
12004). These results give reasonable confidence in using 
local galaxy templates to estimate the infrared luminosi- 
ties of distant IR sources. 

5.2. The total infrared luminosity of the MIPS sources 

at z < 1 

For each library of SEDs previously discussed, we com- 
puted the redshift-dependent relations between Lir and 
the flux density observed at 24 /im. For every source of 
the sample, an estimate of the total IR luminosity was 
thus obtained from each set of templates. These mea- 
sures were weigthed based on the accuracy of their asso- 
ciated library to reproduce the observed relation between 
Lipi and L15 at the derived luminosity (see Fig. 8). They 
were combined to provide the final infrared luminosity 
of the object, and their associated rms was taken as an 
estimate of the corresponding uncertainty. 



10 



E. Le Floc'h et al. 



As noted bv lChapman et all l)2003bj) there might also 
be a bias of MIPS 24 sources toward luminous galax- 
ies with hot dust temperatures compared to the whole 
population of IR-luminous objects at high redshift. This 
could result in additional systematics affecting these lu- 
minosity estimates, and we have not formally considered 
this effect here. We believe however that the uncertain- 
ties due to the SED libraries dominate the systematics 
in our study. 

The IR luminosities are displayed as a function of red- 
shift in Figure 9. Based on the previous arguments we es- 
timate that they are accurate within a factor of ~ 2-3 up 
to z^l. We see that most sources detected below z ~ 0.5 
are only modest infrared emitters (i.e., Lir<1O^^L0), 
and the number of infrared-luminous galaxies found at 
these low redshifts is quite low. At larger distances 
however, LIRGs represent a significant fraction of the 
MIPS galaxy sample in terms of detection number. This 
large population of LIRGs at 0.5 <z< l.O has also been 
observed in various surveys performed at 15 iiia with 
ISOCAM (e.g.. Russel et al., Ji)99; Flores et al. 1999; 
lElbaz et al.M2002|) . We finally note that the most lu- 
minous galaxies (i.e., Ljr > 10^^ Lq) are still pretty rare, 
but there is a clear hint for an increase of the number 
of the brightest LIRGs (5x 10" Lq < Ljr < 10^^ Lq) at 
z ~ 1 , which could point to a large po pulation of ULIRGs 
at even larger distance s (i.e., z > l.fi. lBlain et aT]ll999al: 
iChapman et al.l |2003aj) . These detection rates are ob- 
viously contingent on the comoving volume sampled at 
each redshift though. This general issue will thus be ad- 
dressed in more detail in Sect. 7, where we characterize 
the evolution of the infrared luminosity function based 
on this sample of MIPS sources. 

As a final remark, it is worth mentioning another possi- 
ble caveat related to these infrared luminosity estimates. 
The various template libraries that we used in this 
work are only representative of normal and starburst- 
like galaxies, and do not include the SEDs typical of ac- 
tive galactic nuclei which are significantly flatter in the 
MIR wavelength range. A careful investigation of the na- 
ture of the emission process dominating the 24 /im flux 
in our sub-sample (star-forming activity versus AGN) is 
beyond the scope of this paper, but we note that the 
fraction of MIPS sources showing evidence for the pres- 
ence of an AGN in their optical counterparts is less than 
~ 15% according to the VVDS and COMBO-17 classifi- 
cations. Furthermore, recent synthetic models connect- 
ing the X-ray and infrared SED of AGNs as well as their 
contribution to the cosmic backgrounds in these energy 
bands indicate that the emission arising from pure AGNs 
should be negligible (i.e., < 10%) in high re dshift sources 
detected at 24 //m Ce.g.. iSilva et alJl2004() . In the lo- 
cal Universe moreover, AGNs dominate the mid-IR out- 
put of gala xies o nly for the m ost extreme ULIRGs (e.g., 
iLutz et al.lfl998HTran et al.ll200H) . while the X-ray and 
optical spectral properties of infrared galaxies detected 
with ISO and MIPS at 15 /zm and 24 fim also argue for 
a populatio n dominated by star form ation a t the ^ 85- 
90% level llFadda et all l200llrVa^e schini erani200l 
[Manners e t al.ll2004t lFranceschini~alii2005l : iBell et all 
1^05) . It should therefore be reasonable to assume SEDs 
typical of star-forming galaxies for this current work. 



6. THE OPTICAL COUNTERPARTS OF THE MIPS 
SOURCES: SOME PROPERTIES 

We describe in this section a few basic characteristics of 
the MIPS source optical counterparts to provide a first- 
order answer to very simple questions: are the distant 
infrared luminous galaxies detected by MIPS also lumi- 
nous in the optical or are they highly obscured ? Are 
these dusty starbursts associated with very young sys- 
tems or already evolved galaxies ? A description of their 
morphologies, colors and specific star formation rates is 
provided by BeU et al. (2005). 

The COMBO-17 i?-band magnitudes of the MIPS 
sources are plotted as a function of the fiux density at 
24 /im in Figure 10a. Not surprisingly, there is a clear 
trend for the fainter 24 /im objects to be associated with 
faint optical sources likely located at higher redshifts. 
Such a trend has already been noted among the popula- 
tion of infrared galaxies detected with ISOCAM at 15 /im 
(e.g., Pozzi ct al. 2004). However there is large disper- 
sion in the relation and for most of the sample the op- 
tical brightnesses can vary by more than 5 magnitudes 
when infrared fluxes only change by less than a factor 
of 10. This clearly points to a very wide range of the 
Lm/Lopticai ratios and a broad variety in the nature of 
the MIPS sources (see also the Appendix). 

We also used the absolute magnitude estimates pro- 
vided by COMBO-17 to derive the intrinsic luminosities 
of the MIPS sources at optical wavelengths. As an exam- 
ple we fllustrate in Figure 10b the COMBO-17 absolute 
magnitudes estimated in the B-band filter and corrected 
to our cosmology^ as a function of the total infrared lu- 
minosities derived in the previous section. We observe a 
clear relationship between the two quantities, though the 
dispersion remains relatively large (Icr ~ 0.5 dex). Simi- 
lar relations are also obtained between the total IR lu- 
minosity and the absolute magnitudes derived from the 
other broad-band filters of COMBO-17. They show that 
distant IR luminous sources (at least up to 2;~1) are 
preferentially associated with luminous optical counter - 
parts as in the local Universe fSand ers fc Mira befiggG"). 
This result con firms pre vious studies based on ISO sur- 
veys (jR igopoulo u et al.l 2002: Franc eschini et all l2003t 
iZheng et alj.2004j) . It also suggests that the dust re- 
sponsible for the IR excess in such distant IR-luminous 
objects is likely distributed within small-size regions and 
the corresponding reddening is not sufficient to com- 
pletely obscure the underlying galaxy. This interpreta- 
tion is reinforced by the range of optical colors observed 
among the MIPS 24 /im source population, very simi- 
lar to the range characterizing the local normal galaxies 
(BeU et al. 2005). 

COMBO-17 also provides an estimate of the rest-frame 
galaxy luminosity at A=2800A, which can be used to de- 
rive the "IR excess" in the MIPS sources. IR excess 
is usually defined as the ratio between the IR and UV 
emission and it is plotted in Figure 10c as a function of 
the total IR luminosity. Owing to the tight relation be- 
tween th e star-forming a ctivity and the IR emission of 
galaxies ijKennicuttI 1 1 9981) . this plot can also be read as 
the IR excess as a function of an "IR-equivalent" star 
formation rate (an IR luminosity Ljr = 10^"'^ Lq typically 

^ The COMBO-17 catalog provides absolute quantities assuming 
Ho = 100 km s-l Mpc-i, =0.3 and Qx = 0.7. 
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Fig. 9. — Total infrared luminosities of the MIPS 24 sources identified with a redshift at < z < 1.2 (+ symbols). They were derived 
using various luminosity-dependent SED templates of the literature (see text for more detai ls) and we es timate that they are accurate 
within a factor of 2-3. The median uncertainty of the photometric redshifts is less than 0.05 IWolf et all20d4D . Error bars are not shown 
for clarity. Assuming the calibration from Kcnnicutt (1998), IR luminosities are translated into an "IR-equivalent SFR" reported on the 
right vertical axis. The thick solid line indicates as a function of redshift the infrared luminosity corresponding to an observed 24 fim flux 
of 0.08 mjy (80% completeness limit of our survey). Inset: Corresponding IR luminosity histogram of the sample, with the contribution 
of sources at z < 0.5 (respectively, z < 0.8) indicated by the dark (light) shaded region. 



corresponds to a formation of ~ 17 Mq of stars per year 
following standard calibrations). It illustrates the very 
well-known trend for galaxies characterized by more in- 
tense star-forming activity to be gene rally dustier and 
more luminous at IR wavelengths (e.g..lBuat et al.ll2002l: 
ICardiel et al.ll2003tlPo7zi et al.l l2003': Tlo res et al.ll2004j) . 
We also see that the fraction of UV photons not absorbed 
by dust in the case of the most luminous sources (i.e., 
LjR ^ 10^^'^ L0) is negligible compared to the energy re- 
processed in the IR. 

We finally provide an estimate of the stellar masses of 
the MIPS sources based on the properties of their optical 
counterparts. A tight correlation exists between the rest- 
fram e optical colors of galaxies and t heir mass-to-li ght ra- 
tios ()Bell fc de Jon g 2001; Kauffman ~et al.ll2003|) . This 
relation is uncertain by a factor ^0.4dex due to com- 



bined effects of metallicity, dust, and history of star for- 
mation in individual sources. It is however accurate 
enough for the purpose illustrated here (i.e., getting a 
qualitative estimate of the distribution of IR-luminous 
objects as a function o f mass). F o llowin g the detailed 
analysis and recipes by IB ell et all ()200.1L see their Ap- 
pendix A), we used the y-band absolute luminosities of 
galaxies from COMBO- 17 (converted to our cosmology) 
and transformed these luminosities to stellar masses us- 
ing rest-frame B — V colors. Figure lOd illustrates the 
corresponding histograms for both samples of optically- 
selected sources (thin solid line) and of galaxies detected 
at 24 /im (thick solid line). It is very clear that IR 
galaxies detected by MIPS are on average more mas- 
sive (i.e., Ai'>10^'^ A4q) than the field population of 
sources selected at optical wavelengths. Furthermore, 
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Fig. 10. — a) Observed 24 flux densities as a function of the i?-band magnitudes obtained from COMBO-17. b) Total IR luminosities 
as derived in Sect. 5 versus the absolute B-band magnitudes, c) IR to UV luminosity ratios as a function of the total IR luminosities, d) 
Histogram of the stellar masses (in solar units) for the sample of MIPS sources (thick solid line) compared to the stellar masses of the 
COMBO-17 optically-selected galaxies (thin solid line, scaled down by a factor of 2.2). The dotted-line (respectively, dashed-line) histogram 
corresponds to the stellar mass distribution restricted to 24 /^m sources with Ljj^ > lO'^^ Lq (Ljjj^ > 10^^'^ Lq). Typical uncertainties are 
indicated with error bars in panels b), c) and d). 



the 24 /iiTL objects tend to be associated with even more 
massive galaxies when their IR luminosity (and thus 
their star formation rate) increases. As already observed 
in the case of 15 ^m sources detected by ISO at z < 1 
(|Franceschini et al.ll2003|) . LIRGs with Lir>10"-5Lq 
have most often counterparts with M>1Q^^ M.q. In 
Sect. 8 we will discuss this issue within the more general 
context of the comoving star- forming activity at z ~ I. 

As a side note, we found a group of 91 sources (among 
which 54 have /24pm ^ 83 /^Jy) identified with photomet- 
ric redshifts and totally escaping the various relations 
illustrated on Figs. lOabc. These objects are rather faint 
in the optical (i?>23mag), and their f2ip,m/fR flux ra- 
tios are particularly high and typical of LIRGs/ULIRGs. 
According to our previous results, and given that most 
of them have disturbed merger morphologies c hatacteris- 
tic of high redshift star- forming galaxies fe.g.. lBell et al.l 
1200 5). they are thus likely located at z > 0.5-0.6, which is 
also supported by their morphologies since many of those 
are chatacteristic of high redshift star-forming galax- 
ies. However they have been assigned very low red- 
shifts by COMBO-17, mostly in the apparent overden- 
sity observed at z~0.15 in the CDFS. This would im- 
ply faint infrared luminosities (10^ L0 < Ljr < 10^° L0, 



see Fig. 9) and optical absolute magnitudes much higher 
{Mr > —15) than expected from the correlations. We be- 
lieve that these redshifts have probably been misidenti- 
fied. We decided to exclude the corresponding sources 
from our sample. 

7. INFRARED LUMINOSITY FUNCTIONS 

7.1. Methodology 

We explore in this section the evolution of IR luminos- 
ity functions (LFs) at 0<z<1.2. These LFs were de- 
rived for various r edshift bi ns with the usual 1/V^ax for- 
malis m (ISchmidtlfigea iHuchra fc SargentI 119731 IFeltenI 
11976ft applied to our sub-sample of sources brighter than 
/24/im = 83 /Lt Jy and i? = 24 mag. The selection function 
that was used in this goal was computed as follows. First, 
we considered the probability for a given object to verify 
our selection criterion at 24 /im. This probability equals 1 
for sources brighter than ~ 300 ^ Jy but decreases at 
fainter fluxes due to the incompleteness of the 24 /zm cat- 
alog. The corresponding effect was quantified using the 
detailed simulations described by Paoovich et al. ( 200^ 
see their figure 1). Second, we estimated the probabil- 
ity for a 24 /im source to be associated with an optical 
counterpart at i? < 24 mag. This probability is a redshift- 
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dependent function, and it can be derived with an esti- 
mate of the absolute magnitude corresponding to the op- 
tical selection criterion {R = 24 mag) combined with the 
relation between the infrared and the optical luminosi- 
ties estabhshed in the previous section. At this step, we 
actually assumed that 100% of sources with R< 24 mag 
are detected by COMBO-17^. Finally, we took into ac- 
count the probability for a given optical source to be 
characterized by a redshift. This was done using the es- 
timate of the redshift identification completeness that we 
determined in Sect. 3.3 for sources at i? < 24 mag. 

The computed selection function was subsequently 
used for each single galaxy of the sample to weight the 
differential volume elements that are integrated when cal- 
culating the total accessible volume Knax- Given the 
multi-wavelength flux limits considered here, this inte- 
gration was performed up to the maximum redshift en- 
abling the detection of the object at both 24/L(m and op- 
tical wavelengths (i.e., the lowest value between the red- 
shift where the 24 /im flux would drop below 83 /iJy and 
the redshift where the i?-band magnitude would reach 
24 mag). 

The determination of the luminosity function was per- 
formed for each redshift range independently of the 
sources lying in the other bins. For this reason, the 
Knax approach is more subject to density fluctuations 
than other methods like th e step wise maximum likeli- 
hood re.g.. lEfstathiou et alllTosat IWillmedlTool . Even 
though cosmic variance is not negligible in CDFS, we 
will see however that the uncertainties affecting our con- 
clusions are largely dominated by the errors on the k- 
corrections that are used to translate fluxes into lumi- 
nosities. The Knax technique might also be sensitive to 
the photometric redshift uncertainties, but as we already 
pointed out in Sect. 3.2, the effect should be negligible 
here given the accuracy of the COMBO-17 redshift clas- 
siflcation. 

7.2. Results 

From the luminosity- dependent templates previously 
discussed, we converted the 24 /im observed fluxes into 
monochromatic luminosities at 12, 15, 25 and 60 /im as 
described in Sect. 5.2. The latter were subsequently used 
to derive the associated luminosity functions that we fi- 
nally compared to the local LFs derived from IRAS and 
ISO. In Figure 11 we illustrate the 15/im LF ■i/;i5(L,z) 
estimated for different redshift bins. The corresponding 
data points are reported in Table 1. At this wavelength 
the fc-corrections needed to convert 24 fim fluxes into lu- 
minosities are minimized for most of the sources. The 
majority of sources are indeed located at 0.5<2;<0.8 
where the MIPS 24 /im filter probes the 13.5-16 /im rest- 
frame wavength range. This behavior reduces the de- 
pendence of the luminosity functions on the assumed IR 
templates. We also show in Figure 12 the total IR LFs 
V'ir(L,z) computed from the infrared luminosities shown 
in Fig. 9. Data points are given in Table 2. These LFs 
are obviously more dependent on the choice of IR SEDs 
since the uncertainty in the conversion between flux and 
luminosity is larger ('^0.4dex). They are however easier 

^ This may not be true in the few regions closely surrounding 
very bright stars of the optical image. However this effect is obvi- 
ously negligible over the total field of view considered in this work. 



to interpret in the context of IR galaxy evolution. 

To correct for cosmic variance effects in each redshift 
bin, the luminosity functions were normalized by the ra- 
tio between the _B-band luminosity densities produced 
by galaxies of the blue sequence and measured respec- 
tively in the CDF S and ov er the 3 flelds of COMBO-17 
(iWoff et alJ 120031 12004: Be U et al.ll200l . Even though 
the B-band and the IR selected populations differ in 
their evolution with redshift, they both trace instan- 
taneous star-forming activity (yet with different time 
scales). This connection results in an obvious relation 
between the two that we already highlighted in Sects. 3.4 
& 6 (see Fig. 10b, see also Bell et al. 2005). The charac- 
terization of the i?-band luminosity function over differ- 
ent fields can therefore be used to estimate a first-order 
correction and minimize the cosmic variance affecting our 
24 /im sample. With its redshift peak at z^O.65, the 
CDFS is particularly subject to this effect. 

7.3. Uncertainties 

For each luminosity bin in a given redshift range, un- 
certainties a(Li,Zj) were estimated as the combination 
of (i) poisson noise statistics on the number of sources 

used in the measurement (i.e., the rms error \J^V~^^ y., 

with "k" the source index in the bin) and (ii) uncertain- 
ties in the flux density at 24 /im and the conversion into 
luminosities. The effect of the latter was simulated with 
a Monte Carlo approach. We assigned to each galaxy a 
range of luminosities (monochromatic or integrated) that 
were calculated by taking account of the nominal flux 
and uncertainty at 24 /im as well as the different possible 
SEDs for the conversion. We found this conversion from 
flux to luminosity to be by far the dominant source of un- 
certainty, which explains the larger error bars found in 
the case of i/'ir than for ipi^. Regarding the monochro- 
matic 15 /im LF however, flux uncertainties at 24 /im also 
contribute significantly, especially at the faintest levels. 
For each redshift bin finally, we uniformly added another 
uncertainty (0.15dex) to the whole luminosity function 
based on our estimate of the cosmic variance effect. Be- 
cause it should be negligible, we did not simulate the 
effect related to photometric redshift uncertainties. 

7.4. Parameterization and evolution with redshift 

In both Figures 11 & 12 we have also illustrated the 
luminosity functions determined locally from IRAS and 
ISO. T he 15/im local LF (j)i5{L) was taken from 
iXu et alJ fl998) and Xu (2000). Regarding the total 
IR luminosity function 0ir(L) at z = 0, we show on 
one hand a rece nt result from the r evised IRAS Bright 
Galaxy Sample l|Sanders et al.ll200^ . on the other hand 
an estimate that we derived from the 60 /im local LF 
fSamidcrs et al. 1990; Takcuchi et al. 2003) assuming the 
tight correlation between the 60 iim and the total IR lu- 
minosity of galaxies. Not surprisingly, these two esti- 
mates agree well with each other. A commonly known 
characteristic of these IR luminosity functions is their rel- 
atively large number of sources at the bright end. As a 
result they are better described by a double-exponential 
profile rather than a classical Schechter parameteriza- 
tion. Their usual analytical form is given as follows: 



14 E. Le Floc'h et al. 

TABLE 1 

Luminosity functions derived at 15/iM from the l/Vmax analysis" 



Log [Li5^^ (L0)] $ [# Mpc-» (Log Lu,^m)-^] 



0.3<z<0.45 0.45<z<0.6 0.6<z<0.8 0.8<2<1.0 1.0<z<1.2 

9.0 s.sol^;^ X 10"^ 

9.5 2.93li 30 X 10"'^ 2.62ti i6 x lO^^ 

10.0 5A4+111 X lO-"* l-75lQ7g ^ 2.90l^-29 ^ lO'^ 2.27tJ Jf x lO'^ 2.38ti'^l ^ 1°"^ 

10.5 8-37t6^28 X 10"® 3.781 J;^^ x lO"" 7.72l||^ x lO"* l-40t°;^| x IQ-^ l-SllasI x 10-^ 

11.0 ' ^■'^tl%l X 10-5 3.35l^;^* X 10-5 1.06+°;^? x 10-* 2.54l^;J* x lO-'' 

11.5 ' 8-38l5^29 X 10-*^ l-07io5^79 X 10-5 2.35lf;^* x 10-5 



^assuming a ACDM cosmology with Ho = 70 km s ^ Mpc ^, flm = 0.3 and fix = 0.7. 

TABLE 2 

Total IR luminosity functions derived from the l/Vmax analysis" 



Log [LiR (Lq)] $ [# Mpc-3 (Log Lir)-1] 



0.3<2<0.45 0.45<2<0.6 0.6 < 2 < 0.8 0.8<2<1.0 1.0<z<1.2 

10.0 8-88t"92 X 10"'^ 

10.5 3.52^2:35 X 10-3 3.84l2.5g X 10-=^ 

11.0 0.96tl ll X 10-3 2.29t3;52 x IQ-^ 3.431:2.29 x IQ-^ 2.77lf x 10-^ 

11.5 0.16lg;i3 X 10-3 0.531^^5 x 10-3 1.18t,V " x 10^3 1.78l^;f^ x 10-3 l.lllj',;^* x 10-3 

12.0 6.87+5^15 x 10-5 S.3Ste^29 x 10^5 2.86t'i ll x 10-^ ^■^'^-Iah x 10"" 

12.5 l-7ll?:2| X 10-5 8.38lg^J X 10-6 4.23t3°i^ x 10-5 5.88l4;4" x 10^5 

13.0 l-18lo;88 X 10-5 



^assuming a ACDM cosmology with Hq = 70 km s ^ Mpc ^, Clrn = 0.3 and Qx = 0.7. 
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where dV is the differential element of comoving volume 
and dN(L) the number of sources with a luminosity L 
within dV and per bin of rflogio(i). 



As expected, the comparison between iIj\{L,z) and the 
local luminosity functions reveals a very strong evolution 
of the LFs with lookback time. For each redshift bin, 
iIj\(L,z) was fitted assuming a monotonic evolution of 
the local LF 4'\{L) in both luminosity and density as a 
function of (1+.2): 



i^x{L,z)=g{z)cl>x {L/f{z)) 
with g{z) = (1 + z)°'^ and f{z) = (1 + z)"^ . 



(2) 



We explored the possible range of values for the fitting 
parameters ajy and ai^ using a minimization. This 
was computed from the difference between the observed 
and the fitted luminosity functions over the 0.3 < ^ < 1.0 
redshift range: 



Li 0.i<Zj<l 



i'\,obs. - i'X,model 



(3) 



Because of their lower statistical significance, LFs de- 
termined at z < 0.3 and z>l were ignored in this pro- 
cess. Reduced values were finally transformed into 
likelihood estimates using log(i) = —0.5%^. The corre- 
sponding 1(T, 2(T and 3(t iso-probability contours for -015 
and '0IR a-re displayed in the insets of Figures 11 & 12 
respectively. 

The elliptical shape and the orientation of these con- 
tours in the luminosity and density evolution parameter 

space reflect a well-known degeneracy often encountered 
when fitting high redshift luminosity functions. Given 
the sensitivity of our 24 /im survey indeed, the very faint- 
end portion of ipxi^jz) cannot be constrained at any red- 
shift. The LF can thus be described with a variety of 
scenarios combining different amount of evolution in ao 
and ttL- The constraints that we obtained can be sum- 
marized as follows: 



ipi5{L,z) : 
tpiR{L,z) : 



aL = 2.6 ±0.9 (3cr) 

an = (-1.80 xaL + 6.75) ± 1.1 

ql = 3.15 ± 1.6 (3cr) 

an = (-1.55 x az, + 5.90) ± 1.6 



We will analyze in Sect. 7.4 how this degeneracy can 
be broken by taking account of other independent con- 
straints. Nonetheless we stress that the knee of the lumi- 
nosity functions is well detected throughout the redshift 
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Fig. 11. — The 15 /im luminosity function estimated per redshift bin between z = and z = 1.2 (* symbols) with the l/Vmax formalism. 
The 3cr uncerta inties are indicated with vertical solid lines (see text for more details). Data points are fitted by evolving the local 15 ^tm 
LF (taken from IXi3l.200g and represented as a dotted line in each panel) both in luminosity and density (a^) =2.1, =2.6, solid line). 
Vertical dash-dotted lines denote the luminosity corresponding to an observed 24 /im flux of 83 /ijy (i.e., 80% completeness of the MIPS 
survey) calculated at the lowest redshift considered in each panel. The inset represents the Icr, 2(t and Scr iso-probability contours of the 
likelihood estimated as a function of ajj and with a test. 



range considered in this work (i.e., z<l). This allows 
us to exclude a pure evolution of the LF in density (i.e., 
aL = 0) with a very high confidence. 

Because of the quasi-linearity between L i.»s and Ljr 
(jCharv fc Elbazll200lt Ffakeuchi et alJl2005|) . these con- 
straints on the evolution of "015 and ipm should be in 
principle strictly similar. While they do agree rather well 
with each other, the extension of the contours reveal how- 
ever that the evolution of the 15 LF seems to be bet- 
ter constrained than the evolution of "(Air- As we already 
pointed out, LF uncertainties are largely dominated by 
the errors in the conversion between the 24 /im flux and 
the luminosities, and these errors are in fact minimized 
around 15 /im. Furthermore, we see that ?Ai5 might be 
characterized by a stronger evolution in density and a 
smaller evolution in luminosity compared to V'lR- This 
discrepancy could be a consequence of the non-negligible 
dispersion that exists in the correlations observed be- 
tween fluxes at mid- and far-infrared wavelengths (e.g.. 



HeioullTQM 1x11200(1 IDale et alJ l200lt IDale fc Helol 
200^ Indeed the templates that we used to extrap- 
olate fluxes to luminosities do not really show a pure 
linear relation between L15 and Ljr (see Fig. 8). More- 
over, the constraints on and ul strongly depend on 
the assumed local luminosity functions, and the compar- 
ison between ^15 and 06O (that we used to compute 
clearly reveals a non-linearity between the two. Finally, 
it might be suggested that our scenario based on a single 
evolving population is overly simplistic. Considering dis- 
tinct object_ty£es_cliaiacterizedby_t^^ own evolution 
("Xu"2000f, 'La gache et alll20ml 12001 IPozzi et alJl200l 
Gruppioni et alJl2005t), b i-variate LFs llChapmanet'ld] 
l2003bt iLewis et al.l l2005j) o r luminos itv functions only 
evolving at their bright end l)Charv fc Elbaz 2001) could 
be one way to explore this effect in more detail. 

In spite of this apparent degeneracy, it should be how- 
ever noted that the evolution of the total luminosity 
density fiiR integrated from these infrared LFs is more 
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Fig. 12. — Evolution of the total IR luminosity function (same legend as in Figure 11). The fit has been obtained by evolving the total 
IR local LF (shown as a dotted line in each panel) with ao = 1.0 and oj^ = 3.15. This IR local luminosity functio n was es timated from the 
60 ^im LF of Takeuch i et al. 12003, see also Saunders et al. 1990) and the 60 ^tm/total-IR IRAS correlation Ce.g.. lCharvl; E l baz 200ll). It 
is com pared on the first panel with the total-IR LF derived from the IRAS revised Bright Galaxy Sample (open diamonds. lEanders et alJ 
120031) . The solid grey line represents the total IR luminosity function from the model by Lagache et al. (2004!). In the inset, the dotted 
and dashed lines delimit a region where the parameters ao and lead to an evolution of the star formation activity that is consistent 
with the constraints determined by Hopkins (20(3 see text for details). The shaded region corresponds to the excluded parameter space 
that leads to an overproduction of the 24 fim counts at faint fluxes. This constraint reduces significantly the observed degeneracy. 



tightly constrained than the evolution of the character- 
istic parameters and (j)^ considered separately. In 
each redshift bin, MIPS can indeed detect those sources 
responsible for the bulk of the comoving luminosity den- 
sity, and the additional uncertainty related to the ex- 
trapolation for taking account of the contribution of faint 
objects (i.e., faint-end slope of the LF, see Sect. 7.5) is 
therefore not dominant. Based on the relation that we 
obtained between ol and we find that f2iR evolves 
as (l-|-z)3-9±o-4 at < z < 1 (see also Sect. 8 and Fig. 14). 

Finally, we tested the effect of our redshift identifi- 
cation incompleteness at z>l on the LF estimates by 
computing including the observed luminosity function 
determined at 1 < z < 1.2. For both ipi^ and V-'ir the min- 
imization of in this case leads to a short translation 
of the iso-probability contours toward smaller values of 
ao and larger values of compared to what we previ- 



ously found. This trend is consistent with expectations 
if V'A is underestimated at z>l. Indeed the latter ar- 
tificially yields a fainter evolution in density, the effect 
of which gets compensated in the lowest redshift bins 
thanks to a slight increase in luminosity. As a sanity 
check we finally computed also excluding the LF ob- 
tained at 0.8 < z < 1 where we may also be missing a few 
24 /xm sources fainter than R^2A mag. The results that 
we found remain consistent with those initially obtained. 

7.5. Breaking degeneracies 

iHopkinsI l|2004D combined constraints from source num- 
ber counts at radio wavelengths with estimates of the co- 
moving star formation rate (SFR) density of the Universe 
at high redshift to break the degeneracy that also arises 
when quantifying the evolution of star-forming radio- 
selected galaxies. Following his method we investigate in 
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this section how similar considerations can help in better 
constraining the evolution of IR galaxies. 

Up to z~ 1, the integrated star formation density per 
comoving volume of the Universe p^^fz) is now we ll estab- 
lished within a factor of 2 to 3 (Hopkins 2004) . Given 
the relationship between the total IR emis sion of galaxies 
and their obscured SFR ijKennicuttil 1 998(1 . the evolution 
of V'lR converted into a history of the total dusty star- 
forming activity must be therefore consistent with these 
constraints on Pi,{z). We computed the "IR-equivalent 
SFR" as a function of redshift with different combina- 
tions of a_D and a^. The dispersion in the relation be- 
tween the flux at 24 /im and the total IR luminosity ob- 
viously affects this estimate. Moreover we stress that it 
should only be a lower limit on the true SFR given the 
un- absorbed UV photons produced by young stars and 
not a ccounted for by the IR measurements (e.g., iBelll 
120031 see also the Lir/Luv ratio of the MIPS sources 
in Fig. 10c). Given these caveats and the additional un- 
certainty on Pi,{z) mentioned above, we required our IR- 
SFR determination to lie within 0.5 dex of the averaged 
relation b etwee n logio(/9*(z)) and logio(l+2:) derived by 
iHopkinj 1(20041) . This constraint demarcates a specific 
region in the luminosity and density evolution parame- 
ter space, lying between the dashed and the dotted lines 
shown in the inset of Figure 12. As we can see, the 
evolving parameters au and that we derived from 
minimization in the previous section agree well with 
the known history of star formation up to z 1. However 
we also note that this approach does not really help in 
solving the aforementioned degeneracy. One may need a 
more accurate determination of p-k{z) to progress in this 
direction. Alternatively, a better constraint on the con- 
version between /24^mi Lir and SFR might provide in 
the future a more accurate determination of Pi,{z). 

A much more interesting constraint can be obtained 
from the source number counts. In fact we only con- 
sidered 24 /xm sources brighter than 83 fiJy when build- 
ing the luminosity functions ipx{L,z), but the counts at 
lower fluxes can be used to derive limits on the contribu- 
tion of sources at luminosities fainter than those taken 
to minimize the x^- Using the 24 pra obser vations of the 
"GOODS t est-field" centered on E LAIS-Nl. lGhTrv et a,].l 
((2.004) and iPapovich et alJ l(2004ft constrained the faint 
source density down to f24fj.m ~ 30 pJy and showed that 
the differential counts normalized to the Euclidian slope 
are likely dropping very rapidly below this limit. We 
must therefore ensure that the evolution of the infrared 
LF does not lead to an overproduction of these counts at 
faint fluxes. To check the latter we derived as a function 
of ao and the differential number counts that would 
be produced up to z = 1 by a population of sources de- 
scribed by ip\{L, z) and evolving as given by Equation 2 
(see Fig. 13). We rejected the solutions overproducing 
the total counts obtained from the GOODS test field at 
the faint end. Surprisingly we found that this constraint 
is nearly independent of the library of IR SEDs used to 
describe the galaxy population. This excludes a region 
of the [a£),aL] parameter space that is illustrated by the 
shaded area in the inset of Figure 12. 

This additional constraint reduces significantly the 
number of possibilities to describe the evolution of 
ipYfi^L , z) . The best parameters quantifying this evolu- 



tion at 0<z< 1 are given by cx (1 + z) -« ^ and 

oc (1 -I- z)°'^-o 6 (quoted uncertainty of la) . They 
are summarized in Table 3. In particular they exclude 
any solution favoring a larger evolution in density than 
in luminosity. This trend has already been reported by 
several groups using similar analysis with the ISO and 
SCUBA number counts, and the overproduction of the 
background resulting from a too strong increase of 
is a well-known constraint on the backward evolution- 
ary scenarios of IR galaxies. It has also been noted 
with models assuming pure density evolution in the Press 
Schechter formalism (jMould.200,'^'1 . 




0.01 0.10 1.00 10.00 
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Fig. 13. — Simulated 24 /xm number counts produced up to 
z = 1 by a population tied to the local IR luminosity function and 
evolving in luminosity and density according to various scenarios 
(dotted line: a£i=0.0, 0^=4.0 - dashed line: «£i=0.2, Q!j^=3.25 
- dash-dotted line: 5, ai,=2.75 — triple dot-dashed line: 

a£i=2.5, 0^=2.2). They are compared to the total observed 
24 /im source number counts ('*' symbols, solid line and verti- 
cal error bars) constrained and extrapolated at faint fluxes using 
the "GOODS test field" CCharv ct al. 2004; PaDovich ct al. 2003). 
Note that a too large evolution in density clearly overproduces the 
counts at the faint end. 

Obviously we could also use this method to further 
restrict the possible combinations of a^i and by re- 
jecting scenarios that overproduce the counts at bright 
fluxes (i.e., /24/im ^ 0.3mJy) or do not reproduce the 
cumulative counts discussed in Section 4. For instance 
the model predictions represented by the dashed line in 
Fig. 13 arc consistent with the constraints imposed by 
the faint source counts but severely underestimate the 
true contribution of sources at z < 1. However we found 
that reproducing the bump observed at f^Afim '^0.3mJy 
is very dependent on the choice of SEDs, which could also 
explain why the pre-launch models had over-predicted 
the MIPS 24 /im number counts at these bright fluxes. 
Since our goal is to exclude unphysical solutions without 
being subject to the assumed SEDs, we did not consider 
further this information. 

8. DISCUSSION 

8.1. Evolution and shape of the infrared luminosity 
function 

In the previous section we have constrained the knee 
and the bright end of the infrared luminosity function 
^ir(-L, z) up to z 1. As expected our results indicate a 
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TABLE 3 

Parameterization of the evolution of the total IR luminosity function as a function of redshift" 



Redshift range 


'I'm. 




"IR 


0"IR 


z = (local IRAS luminosity function) 


0.89 X 10^2 


1.77 X 10» 


1.23 


0.72 


0.0<z<0.3 


0.98+° ;;^ X 10-2 


2.77to.o8 X 109 


1.23 


0.72 


0.3<2<0.45 


l.llt°o°l X 10-2 


4.9lto.30 X 10^ 


1.23 


0.72 


0.45<2<0.6 


1.20+027 X 10-2 


6.84+2-35 X 109 


1.23 


0.72 


0.6<2<0.8 


1.29+°-J* X 10-2 


9.7t\ * X 109 


1.23 


0.72 


0.8<2<1.0 




13.8+17 X lO'' 


1.23 


0.72 


1.0<2<1.2 


1.50+°-2* X 10-2 


19.011^6° X 109 


1.23 


0.72 



^assuming a ACDM cosmology with Hq = 70 km s ^ Mpc ^ , Qm =0.3 and Qx 



0.7. 



very strong evolution of this LF with lookback time. We 
find that the space density of galaxies with Ljr > 10^^ Lq 
at z ^1 exceeds by more than 100 their density in the lo- 
cal Universe, which is in fairly good agreement with the 
results from ISO. The quantification of the evolution 
depends however very strongly on the mid- and far-IR 
SEDs used to compute the /c-corrections. This points 
to an urgent need for a more accurate characterization 
of infrared spectral energy distributions of normal and 
luminous galaxies. This goal might be achieved by com- 
bining 24, 70 and 160 fim broad-band MIPS imaging with 
infrared spectroscopy from the IRS spectrograph and the 
MIPS "SED mode" . 

Given the limitations due to the uncertainties of our LF 
estimates, we do not find any evidence for a modification 
of the shape of ipm relative to the local luminosity func- 
tion at least in the luminosity range probed with MIPS. 
Obviously the density of sources fainter than the 80% 
completeness limit of the 24 /im survey is not directly 
constrained and we cannot exclude a slight steepening 
of the LF faint-end slope at high redshifts. However, the 
weight of this LF by luminosity shows that the increase of 
energy produced by faint objects in the case of a steeper 
slope would be mostly driven by sources just below the 
knee of ■(/'ir, the faintest ones having a negligible contri- 
bution in spite of their larger number. The corresponding 
effect is then very similar to the one produced by a too 
large evolution of z) in density. As we saw in the 

previous section it would result in an overproduction of 
the number counts at faint fluxes and violate therefore 
the constraints on the background. While the steepening 
of ipiB. at the faintest end cannot be definitely ruled out, 
it must have marginal significance and will not affect our 
further discussion. 

Furthermore we do not notice any obvious break in 
the shape of the luminosity function 'ijjis,{L, z). Contrary 
to our simple scenario that considers the MIPS detec- 
tions as a single class of objects, these breaks can oc- 
cur when the LF is decomposed into the contribution of 
several populations (e.g., starbursts, AGNs, cold galax- 
ies, ...) evolving independently with the redshift. As an 
example, they can be observed in the model proposed by 
iLagache et ahl feOO.S. 2004). Based on a simple decompo- 
sition of high redshift galaxies into normal non-evolving 
sources and starbursts undergoing strong evolution, this 
model had successfully reproduced previous IR/sub-mm 
observations from ISO, SCUBA and COBE, but it fails 



in explaining the faint-end part of our luminosity func- 
tions because of a predicted break that is not observed 
(see Fig. 12). Spitzer provides therefore new constraints 
on this kind of scenario, and our results suggest in this 
case a smooth transition between populations so that 
continuity is still observed in the total luminosity func- 
tion. This effect could be taken into account consider- 
ing^ multi-varj.a,tc luminosity functions ((Chapman et all 
l2003btlLewis et alJfeOOl . 

8.2. Comparison with results from the previous 
long-wavelength surveys 

Over the last decade cosmological surveys conducted 
with ISO mainly at 15 /im, 90 /im and 170 ^m provided 
direct evidence for the importance of infrared luminous 
sources at 0<2;<1, while submillimeter observations 
with SCUBA revealed a very high density of dusty galax- 
ies with LiR > 10^^ Lq at very high redshifts (i.e., z>2). 
Number counts and redshift distributions that were de- 
rived from these surveys as well as the direct measure- 
ment of the far-infrared background by COBE already 
allowed determination of some constraints on the evo- 
lution of the infrared energy density with redshift (e.g., 
[Blain ct al. 1999b; Franccschini ct al. 2001). 

MIPS is however the first infrared instrument with 
good enough sensitivity to obtain a direct measurement 
of the mid-IR luminosity function up to z ~ 1 (but see 
'Pozzi et all 120041: ISerieant et "all 120041 for LFs at inter- 
mediate redshifts). It might therefore be worth compar- 
ing the constraints obtained from this work with earlier 
results from the literature. First, our 24 /zm survey not 
only excludes very clearly the possibility of a pure evo- 
lution of the IR luminosity function in (f)* , but it also re- 
veals a much stronger evolution of in luminosity than 
in density (i.e., ul ^ cto)- This trend confirms previous 
interpretations of long-wavelength surveys and implies 
that the contribution of dusty luminous galaxies was sig- 
nificantly more important in the past. Our evolution 
constraints (i.e., Lf^(x {l+z)~^-^ and ^^^(x {l+z)~°-^) 
show goo d agreement with previously published analy- 
si s fe.g.. | Blain et alJ Il999albl: iFranceschini et"an 120011: 
Lag ache et al.l I2003t IChapmane^nT20o5 ILewis e'^d] 
fedosj) . However our data do not seem to be consistent 
with a luminosity evolution stro nger than ^ 4.5 as 
in several scenarios examined bv lCharv fc Elba3 l|2001fl 
and iXul l|200Q) . unless one adds a decrease in density 
[ao < 0). 
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Our analysis moreover agrees with the constraints on 
the LFs of star-forming galaxies recently derived from 
the radio sub-mJy source number counts and the history 
of star formation. An evolution slightly lower than our 
estimate (a^ — 2.7±0.6,Q!d = 0.15 ± 0.6) h as been 
obtained with this approach bv iHopkin^ l)2004l) . but it 
is still consistent with our results within the large uncer- 
tainties. The small discrepancy could also originate from 
the likely different levels of AGN contamination charac- 
terizing populations selected at 24 /im and at radio wave- 
lengths. 

While our results on the mid- and total IR luminos- 
ity functions are still consistent with each other, we fi- 
nally note that our constraints on rpi^ seem to require 
a stronger evolution in density and a smaller evolution 
in luminosity compared to V'lR- We have not yet un- 
derstood whether this trend is physically real and re- 
veals the need for a more complex evolutionary scenario 
than the one we have assumed, or whether it is artifi- 
cially produced by the infrared SEDs used to compute 
luminosities. Interestingly, the need for a strong evolu- 
tio n of the mid - lR LF in density has also been suggested 
by iPozzi et all i2004D based o n the ELAIS 15 /im sur- 
vevs fsee also iGruppioni et alJ l2005). A direct compari- 
son between our analysis and their results must however 
be done with caution. Pozzi et al. (2004) decomposed 
their sample into the contribution of non-evolving galax- 
ies that dominate the 15 /im counts at low redshift, and 
distant starbursts responsible for the increase of the IR 
energy density. Given the negligible role of these star- 
bursts at z ~ 0, their evolution must be therefore stronger 
than that of the global sample. This may explain why 
the evolving parameters that they derived are slightly 
larger than ours (i.e., = 3.5t.^'.r;, ap = 3.8t^'n)- Note 
that this approach was also used bv lXu et al.l l)200H) and 
led to a qualitatively similar trend. 

8.3. The evolving contribution of infrared luminous 
galaxies to star formation activity up to z ^ 1 

The constraints on the evolution of ipm^L, z) can be 
used to derive the relative importance of galaxies in a 
given luminosity range and how their contribution to 
the star formation history evolves with redshift. Using 
the fit and the parameterization discussed in the pre- 
vious section we compare in Figure 14 the evolution of 
the IR energy density produced by infrared lumi- 
nous galaxies (i.e., LIRGs -f ULIRGs) and fainter sources 
(i.e., LiR, < 10^^ L0). Uncertainties affecting our results 
are still significant due to the degeneracy previously de- 
scribed. Accordingly the range of possible solutions was 
estimated from the 3cr iso-probability contours showed in 
Figure 12 excluding the combinations of ajj and that 
overproduce the counts at faint fluxes. 

This evolution is also represented in terms of an "IR- 
equiva lent SFR" using the calibration from iKennicutH 
l)1998ft . For comparison we show integrated star forma- 
tion rate densities estimated in various redshift bins and 
taken from the literature (see the compilation bv lHopkinsI 
120041 for references). It should be noted that we have not 
estimated the contribution of the AGN IR emission to 
Following the arguments discussed in Sect. 5.2, we 
believe that such contribution results in a 10-15% over- 
estimate in the true SFR. On the other hand this IR- 
SFR estimate does not take into account the contribu- 



tion of the unabsorbed UV light produced by the young 
stars. As a result we are likely underestimating the total 
star formation rate density by a factor ranging between 
^ 20-30% at 2 ~ 1 (where dusty galaxies dominate the 
SFR) and ~ 60-70% at z ~ (where the star-forming 
activity occurs within fainter sources with low extinc- 
tion) . This effect can be seen by considering, in addition 
to the IR emission, the evolution of the UV luminosity 
(uncorrected for dust extinction) with redshift. Using 
GALEX data. lSchiminovich et all l)2005j) found that the 
energy density measured at 1 500A evolves as (1-l-z)"'^-^ 
at < z < 1 , which is also consisten t with previous UV 
measurements fe.g.. lLillv et al.ll996l) . An estimate of the 
total star formation rate density (represented as a dotted 
line in Fig. 14) can thus be obtained by adding the con- 
tribution of the equivalent UV-uncorrected SFR (dashed 
line in Fig. 14) to our best fit of the evolution of the IR 
star formation estimate (converted from the energy den- 
sity iliR evolving as (1-l-z)'^'^). This clearly shows that 
the IR emission provides a good approximation of the to- 
tal SFR density at z ~ 1 given the uncertainties affecting 
its current measurements (i.e., conversion between flux 
and luminosity, cosmic variance). 

One can see that the global evolution of the luminosity 
function leads to a source density increase for both popu- 
lations of IR-luminous and low luminosity galaxies. How- 
ever the larger increase of L^^ compared to results in 
a much more rapid evolution of sources characterized by 
the highest bolometric luminosities. The contribution of 
IR luminous galaxies, though negligible in the local Uni- 
verse, becomes comparable to that of normal starbursts 
around z~0.7 and they dominate beyond. We also see 
that such IR luminous sources are mostly dominated by 
LIRG-type objects up to z 1 and that ULIRGs still 
have a modest impact (~ 10%) at this redshift. This is 
once again in good agreement with the ISOCAM surveys 
fe.g.. lAussel et al.i ri999). However these ULIRGs have 
undergone the fastest evolution in the last 8 Gyrs and 
our results suggest that their contribution is likely still 
rising at z > 1 . As already revealed by the submillime- 
ter surveys fsee lBlain et al.ll2002l for a review), they may 
therefore be responsible for a very significant fraction of 
the star-forming activity of the Universe at z ~2-3 (see 
also Cowie et al. 2004). 

How can this strong evolution observed in the in- 
frared be understood in the general context of the 
growth of structures ? There is now increasing ev- 
idence that the most massive galaxies seem to have 
formed their stars early in cosmic history and that their 
contribution to the comoving SFR density of the Uni- 
verse was s ignificantly larger a t higher redshifts (e.g. , 
Cowie et al.ll]^9l iCimatti et al.ir2 004: Da ddi et alJl2004t 
Juneau et alJl2005^ . Similarly it has been shown that a 
significant fraction of galaxies with masses Ai > 10^" TW© 
at z~ 0.7-1.0 are experiencing a violent episode of star 
formation, while nearly all equally-massive sources at the 
present day are fo und in a very quie scent r node (e.g.. 
[Flores et all 119991: i Frances chini et all 120031 iBell et alJ 
12005: Ha mmer et al . 2005). In parallel to this downsizing 
effect characteristic of the stellar mass assembly history, 
the evolution of the luminosity function at infrared and 
submillimeter wavelengths reveals a transfer of the star 
formation from the most luminous high redshift sources 
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Fig. 14. — Evolution of the comoving IR energy density up to 2 = 1 (green-filled region) and the respective contributions from low 
luminosity galaxies (i.e., Lm < 10^^ Lq, blue-filled area), "infrared luminous" sources (i.e., Ljr > 10^^ Lq, orange-filled region) and ULIRGs 
(i.e., LiR > 10^^ L0, red-filled region). The lower and upper curves delimiting these regions result from the degeneracy in the evolution of 
i/'iR (see Sect. 7 for more details). The solid line evolves as and represents the best fit of the total IR luminosity density at ^ 2 ^ 1. 

Estimates are translated into an "IR-equivalent SFR" density given on the right vertical axis, where an absolute additional uncertainty 
of ~0.3dex should be added to reflect the dispersion in the conversion between luminosities and SFR. Note that the percentage of the 
contribution from each population is likely independent of this conversion. The dashed line corresponds to the SFR measured from the UV 
luminosity not corrected from dust extinction. The dotted line represents the best estimate of the total star formation rate density as the 
sum of this uncorrected UV contribution and the best fit of the IR-SFR (solid line). At 2 ~ 1 IR luminous galaxies represent 70±15% of the 
comoving IR energy density and dominate the star formation activity. Open diamonds, vertical and horizontal bars represent integrated 
star formation rate densities and their unc ertainties estimated within various redshift bins and taken from the liter ature I Connolly et aU 
1997^ Cli-esse & Maddox"1998; 'li-ever et al.''l998; Florcs et al.ll l99a ICowie et alj[T99gt IHaarsma et al.ll200d IMachalski &: Godlowski..200oF 
S^livan et al. 2001; Condon et al. 2002; Sadler et al. 2002; Serje ant et all2002tlTresse et al.l2002tlWilson et al.i2OO2tlPere2-Gon2a.le2 et^ 
2003tiPo2zi et al...2004; see the compilation bv ■Hopkinsi:i,2004) . 



to more modest local starbursts. If we consider that 
LIRGs and ULIRGs are preferentially associated with 
more massive systems (see Fig. lOd), this result might 
be closely related to the mass growth picture described 
above. However the connection between the two is likely 
more subtle as we found that the fraction of IR lumi- 
nous sources with A4 < 10^° is not negligible either. 
The relation between IR luminosities and masses has in- 
deed a large dispersion and the dominant contribution 
of LIRGs/ULIRGs to the star formation at z^l does 
not necessarily imply that the SFR at these redshifts is 
only locked in ma ssive systems. A sim ilar conclusion has 
been proposed bv lJuneau et al.l <|2005j) . who showed that 



the star formation density at z 1 was dominated by 
systems with < 10^°'^ TW©. Much larger samples of 
sources will be required to address this issue in more de- 
tail, especially by de-projecting the comoving SFR den- 
sity into the IR-luminosity/mass/redshift 3-dimensional 
space. This goal could be achieved in the very near future 
by combining all the MIPS cosmological surveys that are 
currently being carried out with Spitzer. 

9. CONCLUSION 

We have analyzed a sample of MIPS/ Spitzer 24 /im 
sources detected in the Chandra Deep Field South us- 
ing ancillary optical data from the literature (i.e., mag- 
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nitudes, spectroscopic and photometric redshifts). Our 
results can be summarized as follows: 

• For 24 /im sources brighter than ~ 80 /iJy, the se- 
lection criterion i? < 24 mag provides a complete 
sample of optical counterpart identifications up to 
z-0.8. 

• About 55-60% of the 24 fj,m objects brighter than 
~ 80 /iJy are located at ^ < 1, which points to a sig- 
nificant fraction of the MIPS 24 /im sources being 
luminous and ultra-luminous infrared galaxies at 
even higher redshifts. 

• The 24 source population at 0.5<z<1.0 is 
dominated by LlRG-type objects and slightly 
fainter sources (i.e., Lir < 10^^ Lir). ULlRGs are 
found to be rare at these redshifts. 

• The conversion between fluxes and luminosities de- 
pends strongly on the assumed infrared SEDs. This 
points to a crucial need for more accurate deter- 
minations of IR templates, a goal which could be 
achieved with MIPS and IRS. 

• As in the local Universe, infrared luminous sources 
up to z ~ 1 are also luminous at optical wavelengths 
and they tend to be more massive than the bulk of 
optically-selected distant galaxies. 

• The comoving energy density measured in the in- 
frared evolves as (l-t-z)^-^^" '' at < z < 1. In con- 
trast, the luminosity density in the UV only evolves 
as (l-t-z)"^^'^ over the similar redshift range. This 
points to a more important extinction by dust re- 
processing light in the IR at high redshift. 

• The infrared-selected sources at < z < 1 have 
undergone strong evolution characterized by 



21 

, o T + 0.7 , , „ 7+0.2 

oc (l+z)^-2-o 2 and (/.^p, oc (l+z)°-^-o e. 
At z^l, infrared luminous galaxies (i.e., 
Lir^IO^^Lq) appear to be responsible for 
70±15% of the comoving IR energy density. They 
dominate the star-forming activity beyond z ~ 0.7. 
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APPENDIX 

SELECTION EFFECTS IN THE SAMPLE OF 24 /iM SOURCES 

In this work we imposed a selection criterion i? < 24 mag when considering the 24 /im sources identified with a 
photometric redshift from the COMBO-17 survey. As we showed in Sect. 3, this selection obviously results in an 
incompleteness of the MlPS-selected population at z> 0.8-1.0. We provide in this Appendix a few more details 
characterizing the possible nature of the sources that we may have missed in the estimate of the LFs at high redshift. 

In Figure 15 we show the relation between the R — 24/im color and the i?-band magnitude of the MIPS sources 
identified with an optical counterpart from COMBO-17 (top panel). The 24/im Vega magnitudes were calculated 
assuming a Zero Point of 7.3 Jy as defined in the Spitzer Observing Manual}'^ As pointed out in Sect. 6, the apparently 
well-defined correlation that is observed is a natural consequence of the small range of 24 /im fluxes (typically a factor 
of ~ 10 over our sample) compared to the much larger range (~ 5 magnitudes) covered by the optical counterparts 
of the MIPS sources (see also Fig. 8a). If we only consider 24/im objects brighter than the 80% completeness of our 
survey (i.e., /24^m > 83/iJy), we see that the selection at i?= 24mag defines a complete sample up to _R — 24/iTO ~ 11.6. 
We miss therefore the reddest sources of the MIPS catalog (shaded region of Fig. 1). 

To get a hint into the possible nature of these targets, we also plot in figure 15 the total infrared luminosities (as 
derived in Sect. 5) versus the R — 2Aiim colors for the MIPS sources identified with a spectroscopic or a photometric 
redshift (bottom panel). Different symbols are used to highlight sources as a function of the i?-band magnitude. As 
we see, 24 /xm objects redder than R — 24/im ~ 11.6 are not obviously the most extreme in terms of dust emission and 
their luminosity can vary between ^SxlO^^L© and a few 10^^ L© in the infrared. Even though sources fainter than 
i?~ 24mag may follow a different trend, this dispersion of Lir with both the color and the i?-band magnitude suggests 
that these sources span a wide range of IR luminosities. Taking into account the evolution of the sensivity limit with 
redshift in our 24 /tm survey (see Fig. 7), they should thus be located in a wide range of redshifts at z > 1. 
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Fig. 15. — Top: The i?-24 /im color as a function of the i?-band magnitude for the MIPS sources identified with an optical counterpart 
in COMBO-17. Sources above the solid line are brighter than the 80% completeness limit of the MIPS 24 /im survey (i.e., /24^m > 83 /^Jy). 
Given the selection criterion i? < 24 mag (dashed line) , sources redder than i?-24/jm~11.6inthe shaded region are therefore missed when 
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